import numpy as np


#计算时域特征
def tdi(data):
    #  有量纲时域特征指标
    max = max(data);                      #最大值
    min = min(data);                      #最小值
    pkpk = max(data)-min(data);           #峰峰值
    mean = np.mean(data);                    #均值
    abs = np.mean(abs(data));                #绝对平均值
    rms = np.rms(data);                      #均方根值
    var = np.var(data,1);                    #方差
    std = np.std(data,1);                    #标准差
    xfg = np.mean(np.sqrt(abs(data)))^2;        #方根幅值
    #   无量纲时域特征指标
    Kurtosisfactor = kurtosis(data,1);           #峭度因子
    Skewnessfactor = skewness(data,1);           #偏度因子
    Peakfactor = max/rms;               #峰值因子
    Marginfactor = max/xfg;             #裕度因子
    Skewness = Skewnessfactor*std^3;   #偏度
    Kurtosis = Kurtosisfactor*std^4;   #峭度
    Pulsefactor = max/abs;              #脉冲因子
    Wavefactor = rms/abs;               #波形因子
    return

def kurtosis(data,type):

    return

def skewness(data,type):

    return